Two phase mixed convection Al2O3–water nanofluid flow in an annulus

Heat transfer enhancement of a mixed convection laminar Al₂O₃–water nanofluid flow in an annulus with constant heat flux boundary condition has been studied employing two phase mixture model and effective expressions of nanofluid properties. The fluid flow properties are assumed constant except for the density in the body force, which varies linearly with the temperature (Boussinesq’s hypothesis), thus the fluid flow characteristics are affected by the buoyancy force. The Brownian motions of nanoparticles have been considered to determine the effective thermal conductivity and the effective dynamic viscosity of Al₂O₃–water nanofluid, which depend on temperature. Three-dimensional Navier–Stokes, energy and volume fraction equations have been discretized using the finite volume method while the SIMPELC algorithm has been introduced to couple the velocity–pressure. Numerical simulations have been presented for the nanoparticles volume fraction (?) between 0 and 0.05 and different values of the Grashof and Reynolds numbers. The calculated results show that at a given Re and Gr, increasing nanoparticles volume fraction increases the Nusselt number at the inner and outer walls while it does not have any significant effect on the friction factor. Both the Nusselt number and the friction coefficient at the inner wall are more than their corresponding values at the outer wall.     

基于汽车与行人碰撞载荷特点的下肢长骨建模

基于THUMS(total human model for safety)下肢长骨有限元模型, 在材料和单元属性等方面进行了改进. 在详细分 析行人下肢长骨载荷特点的基础上, 采用多种不同载荷工况下较新的生物力学实 验数据, 对长骨模型进行了前-后和外-内加载 方向的准静态三点弯曲验证, 同时对近心端1/3处、中部和远心端1/3处加载的动态三点弯曲 验证. 验证结果表明, 该模型具有较好的生物逼真度, 能够准确地模拟行人下肢长骨的骨折 及碰撞响应, 可用于后续行人下肢模型的开发, 并为行人下肢损伤机理和安全防护研究提供 准确高效的研究手段.     

Numerical simulations of counter-current two-phase flow experiments in a PWR hot leg model using an interfacial area density model

In order to improve the understanding of counter-current two-phase flows and to validate new physical models, CFD simulations of 1/3rd scale model of the hot leg of a German Konvoi PWR with rectangular cross section was performed. Selected counter-current flow limitation (CCFL) experiments at the Helmholtz–Zentrum Dresden–Rossendorf (HZDR) were calculated with ANSYS CFX 12.1 using the multi-fluid Euler–Euler modeling approach. The transient calculations were carried out using a gas/liquid inhomogeneous multiphase flow model coupled with a k-ω turbulence model for each phase. In the simulation, the surface drag was approached by a new correlation inside the Algebraic Interfacial Area Density (AIAD) model. The AIAD model allows the detection of the morphological form of the two phase flow and the corresponding switching via a blending function of each correlation from one object pair to another. As a result this model can distinguish between bubbles, droplets and the free surface using the local liquid phase volume fraction value. A comparison with the high-speed video observations shows a good qualitative agreement. The results indicated that quantitative agreement of the CCFL characteristics between calculation and experimental data was obtained. The goal is to provide an easy usable AIAD framework for all Code users, with the possibility of the implementation of their own correlations.     

Flow regime transition criteria for two-phase flow in a vertical annulus

In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug–bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been checked against air–water as well as boiling flow regime maps. In all the cases, the new developed model shows better predicting capabilities than the existing correlations most used in literature.     

A new wall-damping function for large eddy simulation employing Kolmogorov velocity scale

A new wall-damping function, based on the Kolmogorov velocity scale, for large eddy simulation (LES) is proposed, which accounts for the near-wall effect. To calculate the Kolmogorov velocity scale, u_ε, the dissipation rate of turbulent energy, ε, is needed. In LES, however, the dissipation rate is generally not solved, unlike in the Reynolds averaged Navier-Stokes (RANS) simulations, e.g., k-ε models. Although, in some previous studies, the dissipation rate of the subgrid-scale (SGS) turbulent energy, ε_SGS, is used instead of ε in calculating the Kolmogorov velocity scale, the scale obtained using such a method overly depends on the grid resolution employed and is generally inappropriate. Accordingly, the wall-damping function using the incorrect velocity scale also depends on the grid resolution and gives an inadequate wall effect. This is because ε_SGS contains only the components in the scale smaller than the grid-filter width, which obviously varies with the grid resolution employed. In this study, to overcome this problem, we propose a method for estimating the Kolmogorov velocity scale with a technique of conversion in LES, and the estimated one is utilized in the wall-damping function. The revised wall-damping function for LES is tested in channel flows and a backward-facing step flow. The results show that it yields a proper near-wall effect in all test cases which cover a wide range of grid resolution and Reynolds numbers. It is also shown that all three kinds of SGS models incorporating the present wall-damping function provide good predictions, and it is effective both in one-equation and 0-equation SGS models. These results suggest that the use of the proposed wall-damping function is a refined and versatile near-wall treatment in LES with various kinds of SGS models.     

Effect of thermophoresis and other parameters on the particle deposition on a tilted surface

In this study a modified version of v2-f turbulence model (φ-α), is applied to simulate a non-isothermal air-flow. The φ-α model and a two-phase Eulerian approach complement each other to predict the rate of particle deposition on a tilted surface. The φ-α model can accurately calculate the normal fluctuations, which mainly represent the non-isotropic nature of turbulence regime near the wall. The Eulerian model was modified considering the most important mechanism in the particle deposition rate when compared to the experimental data. The model performance is examined by comparing the rate of particle deposition on a vertical surface with the experimental data in a turbulent channel flow available in the literature. The effects of lift force, turbophoretic force, thermophoreric force, electrostatic force, gravitational force and Brownian/turbulent diffusion were examined on the particle deposition rate. The results show that, using the φ-α model predicts the rate of deposition with reasonable accuracy. The results of modified particle model are in good agreement with the experimental data. This study highlights the paramount effect of thermophoretic force on the particle deposition rate and clearly shows that when the temperature difference exceeds a certain limit, the electrostatic force has insignificant effect on the particle deposition rate. Furthermore, it is indicated that even at small temperature differences, the effect of tilt angle on the particle deposition rate for intermediate-size particles is negligible.     

The issues of the uniqueness and the stability of the homogeneous response in uniaxial tests with gradient damage models

We consider a wide class of gradient damage models which are characterized by two constitutive functions after a normalization of the scalar damage parameter. The evolution problem is formulated following a variational approach based on the principles of irreversibility, stability and energy balance. Applied to a monotonically increasing traction test of a one-dimensional bar, we consider the homogeneous response where both the strain and the damage fields are uniform in space. In the case of a softening behavior, we show that the homogeneous state of the bar at a given time is stable provided that the length of the bar is less than a state dependent critical value and unstable otherwise. However, we also show that bifurcations can appear even if the homogeneous state is stable. All these results are obtained in a closed form. Finally, we propose a practical method to identify the two constitutive functions. This method is based on the measure of the homogeneous response in a situation where this response is stable without possibility of bifurcation, and on a procedure which gives the opportunity to detect its loss of stability. All the theoretical analyses are illustrated by examples.     

On the predictive capabilities of the shear modified Gurson and the modified Mohr-Coulomb fracture models over a wide range of stress triaxialities and Lode angles

The predictive capabilities of the shear-modified Gurson model [Nielsen and Tvergaard, Eng. Fract. Mech. 77, 2010] and the Modified Mohr-Coulomb (MMC) fracture model [Bai and Wierzbicki, Int. J. Fract. 161, 2010] are evaluated. Both phenomenological fracture models are physics-inspired and take the effect of the first and third stress tensor invariants into account in predicting the onset of ductile fracture. The MMC model is based on the assumption that the initiation of fracture is determined by a critical stress state, while the shear-modified Gurson model assumes void growth as the governing mechanism. Fracture experiments on TRIP-assisted steel sheets covering a wide range of stress states (from shear to equibiaxial tension) are used to calibrate and validate these models. The model accuracy is quantified based on the predictions of the displacement to fracture for experiments which have not been used for calibration. It is found that the MMC model predictions agree well with all experiments (less than 4% error), while less accurate predictions are observed for the shear-modified Gurson model. A comparison of plots of the strain to fracture as a function of the stress triaxiality and the normalized third invariant reveals significant differences between the two models except within the vicinity of stress states that have been used for calibration.     

Micromechanical modelling of remanent properties of morphotropic PZT

The present paper investigates the capability of micromechanical material models to predict the ferroelectric behaviour of morphotropic PZT ceramics in a rate-independent approximation based on realistic microscopic material parameters. Starting point is a three-dimensional tetragonal model, which builds on the model of Pathak and McMeeking [2008. Three-dimensional finite element simulations of ferroelectric polycrystals under electrical and mechanical loading. Journal of the Mechanics and Physics of Solids 56, 663-683]. Volume fractions of the crystallographic variants represent the domain structure inside the grains. Interactions between the grains are taken into account by means of a representative volume element of the grain compound. A simplified set of realistic microscopic material parameters of the lattice in terms of Young's modulus, Poisson's ratio, dielectric constant, and spontaneous strain and polarisation is derived from experimental data and theoretical results given in the literature. The simulation of the macroscopic remanent polarisation and strain response due to two load cases shows explicitly that the tetragonal model is not capable to reproduce the behaviour of morphotropic PZT. Therefore, the model is extended by the rhombohedral phase, allowing a mixture of both phases with varying quantities inside the grains. A comparison of our results with experimental data shows a remarkably good agreement, revealing the capability of the extended model.     

The interlayer shear effect on graphene multilayer resonators

Graphene nanostrips with single or few layers can be used as bending resonators with extremely high sensitivity to environmental changes. In this paper we report molecular dynamics (MD) simulation results on the fundamental and secondary resonant frequencies f of cantilever graphene nanostrips with different layer number n and different nanostrip length L. The results deviate significantly from the prediction of not only the Euler-Bernoulli beam theory (f∝nL⁻²), but also the Timoshenko's model. Since graphene nanostrips have extremely high intralayer Young's modulus and ultralow interlayer shear modulus, we propose a multibeam shear model (MBSM) that neglects the intralayer stretch but accounts for the interlayer shear. The MBSM prediction of the fundamental and secondary resonant frequencies f can be well expressed in the form f−f_mono∝[(n-1)/n]^bL^−2(1−b), where f_mono denotes the corresponding resonant frequency as the layer number is 1, with b=0.61 and 0.77 for the fundamental and secondary resonant modes. Without any additional parameters fitting, the prediction from MBSM agrees excellently with the MD simulation results. The model is thus of importance for designing multilayer graphene nanostrips based applications, such as resonators, sensors and actuators, where interlayer shear has apparent impacts on the mechanical deformation, vibration and energy dissipation processes therein.